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Superconducting quantum circuits, fabricated with multiple layers, are proposed to implement 
perfect quantum state transfer between nodes of a hypercube network. For tunable devices such as 
the phase qubit, each node can transmit quantum information to any other node at a constant rate 
independent of the distance between qubits. The physical limits of quantum state transfer in this 
network are theoretically analyzed, including the effects of disorder, decoherence, and higher-order 
couplings. 
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I. INTRODUCTION 

Quantum information processing requires a large num- 
ber of highly interconnected qubits. However, most 
leading experimental candidates suffer from limited con- 
nectivity. Typical designs involve low-dimensional net- 
works with nearest-neighbor interactions, such as spins 
of electrons^ or nuclei^ in solids, or atoms in optical 
lattices.— Other designs have networks of a limited num- 
ber of qubits interacting through a common mode, such 
as ion traps. 4 In all these cases, the underlying spa- 
tial arrangement leads to a minimum time for infor- 
mation to be routed between the most distant qubits, 
a challenge to scalability. Moving quantum informa- 
tion through the computer by swap gates, or moving 
the qubits themselves, are plausible but potentially slow, 
error-ridden processes. More sophisticated proposals in- 
clude measurement-based teleportation protocols^ and 
coupling of distant qubits by photons^ both pose addi- 
tional experimental challenges. Efficient quantum rout- 
ing remains a significant design problem for quantum- 
computer architectures. 

Superconducting circuits are a remarkable exception, 
as their couplings are controlled by the fabrication 
of superconducting wires. This wiring can be com- 
plex, involving multiple layers (using interconnects and 
crossovers)^ and thus is capable of three- or higher- 
dimensional topologies. From a quantum information 
perspective, this resource in connectivity can be exploited 
in novel architectures^ and can even be integrated with 
other physical qubits such as atoms^ ions^ or polar 
molecules^ From a physics perspective, studying these 
new artificial solids opens up a number of opportunities. 

The simplest such study is quantum transport, partic- 
ularly the coherent transfer of one qubit state to another. 
As shown by Bose, state transfer in a one-dimensional 
chain, using time evolution of a fixed Hamiltonian, can 
efficiently distribute entanglement over large distances 12 
without photons or feedback. Christandl et al. studied 
networks implementing perfect state transfer— and Feder 
has recently shown how to generate an infinite hierarchy 
of such networks^ One such network is the hypercube, 
well-known in classical computer design^ and even the 
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FIG. 1: (a) Hypercube network of dimension d = 4, with 
corners labeled a and b. (b) Each node is implemented as 
a current biased Josephson junction, while each connection 
requires a coupling capacitor. Typical circuit parameters in- 
clude C x = 6 pF, C c = 60 fF, I c = 21/uA, and w /27T = 6 
GHz. (c) Crossing lines in the network are implemented by 
superconducting wires in a crossover configuration, using ei- 
ther insulating material or vacuum gaps between the layers. 
For sufficiently large gaps (100 nm - 1/nn) the capacitance 
between the wires can be made negligible, (d) A quadratic 
frequency shift of each row of qubits in the hypercube is re- 
quired for high-fidelity state transfer (see text). 



subject of investigation by Feynman one summer years 
ago^ 

We propose to use multi-layered superconducting cir- 
cuits of capacitively-coupled phase qubits 1 -^ to implement 
hypercube quantum state transfer. Coherent manipula- 
tions of single phase qubits using Rabi oscillations have 
been successful in a number of recent experiments, 18 
while coupled systems have also been studied^ includ- 
ing full quantum state tomography of entanglement 31 In 
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addition to direct coupling of two qubits, superconduct- 
ing circuits allow for coupling through a third resonant 
circuit. This was first demonstrated in spectroscopy ex- 
periments with phase qubits coupled to a lumped-circuit 
resonator— and more recently in time-domain experi- 
ments involving the active emission and absorption of a 
photon in a transmission line coupling two phase qubits 22 
and through dispersive interactions of charge qubits with 
an intermediate cavity^ 

While we focus on phase qubits coupled by capacitors 
in the following, our results can be readily extended to 
these and other superconducting qubit designs. Previ- 
ous theoretical studies of one-dimensional quantum state 
transfer have been performed for both charge^ and flux 25 
qubits. More generally, the fundamental importance of 
novel quantum transport in superconducting qubit arrays 
was discussed by Levitov et o/.— Central to our multi- 
dimensional circuit design is the use of long-range cou- 
plings and multi-layered devices — these properties have 
also been used in other recent designs for superconduct- 
ing quantum computers^ 

In the following sections, we first describe the design of 
a phase qubit hypercube. Then, by using the tunability 
of the phase qubit, we show how high fidelity state trans- 
fer can be performed between any two of the 2 d qubits in 
the d-dimensional hypercube network. Furthermore, we 
show how to correct errors due to the higher-order cou- 
pling terms naturally present in the multi-qubit Hamilto- 
nian. Finally, we analyze the effects of decoherence and 
disordered couplings to show that this approach can be 
demonstrated using existing technology. 

II. PHASE QUBIT HYPERCUBE 

We consider a network of capacitively-coupled phase 
qubits— (see Fig. [1]), modeled as current-biased Joseph- 
son junctions described by the following Hamiltonian 

3 

(i) 

where junction j has a critical current I c j and adjustable 
bias current Ij, and where <&o = h/(2e) is the flux quan- 
tum. The dynamical variables are the phase differences of 
each junction jj and their corresponding conjugate mo- 
menta pj , the latter coupled by (the inverse of) the capac- 
itance matrix Cjk- In the ideal case (to be relaxed in Sec. 
V), all junctions have the same critical current I C j = I c , 
intrinsic capacitance C x , and are coupled by identical ca- 
pacitors C c . The capacitance matrix then takes the form 
Cjk — C(5jk — (Ajk), where C = C x + dC c is the total ca- 
pacitance of each junction, ( = C c /C = C C /(C X + dC c ) 
is a dimcnsionless coupling parameter, and Ajk is the 
adjacency matrix for the d-dimensional hypercube (see 
below) . 

By expanding the Hamiltonian in terms of the lowest 
two eigenstates of each uncoupled phase qubit, and uti- 



lizing a rotating wave approximation, we find 

n * -\ E fojZj + \y. m * + ( 2 ) 

3 3<k 

where A, Y, and Z are the Pauli matrices for each qubit, 
u 3 « (27r/ c (C- 1 ) w /$ ) 1/2 (l - and the cou- 

pling matrix is 

n Jk = \^j-ur k ((A + ( 2 a 2 + c 3 a 3 + • • • ). k . (3) 

The adjacency matrix can be written most efficiently 
by first labeling each vertex by a binary number x with 
d bits, i.e. x = x\X2 ■ • • x& with each Xj equal to or 
1. A hypercube will result if we connect vertices whose 
labels differ only by one bit. Thus, the corner labeled 
x = 000. ..0 will have d neighbors, each with a single 1 in 
its label (e.g. y = 100... 0). The elements A X:V of the ad- 
jacency matrix are unity whenever the labels correspond 
to connected vertices, and are zero otherwise. 

Using the binary labeling one can introduce a tensor- 
product structure on the adjacency matrix. Letting r be 
the Pauli matrix r — a x operating in this vertex space, 
one has 

A = 5>> (4) 

3=1 

where each term Tj in the sum corresponds to the ten- 
sor product of r acting on the binary index j with the 
identity matrix for every other index: 

( T j)x,y — T Xj,yj Y\ 3&k,Vk- 

For example, the adjacency matrix of the cube (d = 3) 
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An experimental implementation of this design re- 
quires a fabrication process using multiple layers. The 
hypercube topology can be accomplished with only pair- 
wise crossing of the coupling wires. These wires can be 
fabricated in a cross-over configuration, shown in Fig. 
1(c), in which insulating material separates the two su- 
perconducting layers. By making this separation large 
enough (~ 1/^m), cross-coupling between these wires can 
be reduced significantly. Multi-layered superconducting 
circuits with six metallic layers have been reported^ 

In addition, coupling a large number of qubits will re- 
quire wires spanning multiple qubits. A long wire will 
have a non-negligible inductance L, leading to a reso- 
nant mode with frequency u>lc = 1 / \/LC c which could 
perturb the state transfer dynamics. However, long- 
range coupling of this sort has already been tested^ with 
lengths approaching one millimeter, large compared to 
the typical dimensions of the qubits (~ 10/zm). The cor- 
responding inductance in this circuit was found to be 
L ~ 2 nH, which combined with a coupling capacitance 
C c = 30 fF leads to uj lc /2tt > 20 GHz, much higher 
than the qubit frequency (ujq/2tt ~ 6 GHz). This is sig- 
nificantly off-resonant and thus the resonant mode will 
remain in its ground state — although the nature of the 
coupling will be slightly modified.— 

Finally, the actual fabrication of a hypercube circuit 
appears quite complex. However, the number of wires 
in this design is actually quite modest. The hypercube 
circuit can couple any two qubits in the same amount 
of time — this is shown in the next section. This prop- 
erty is shared by a completely connected set of N qubits. 
However, since a <i-dimensional hypercube has d2 d ~ 1 
edges^ the hypercube requires significantly less con- 
nections. Specifically, the completely connected circuit 
would require N(N - l)/2 wires (~ 10 6 for N = 1024), 
while the hypercube would require only (iV/2) log 2 AT 
wires (~ 5000 for N = 1024). We conclude that this is 
a promising theoretical design that can be implemented, 
at least for modest d, using existing technology. 



III. PERFECT STATE TRANSFER 

Quantum state transfer involves the preparation of an 
initial state = (a|0) +/?|l)) o ®|0) rest ® |0) 6 , contain- 
ing quantum information in the qubit a, and evolution 
by the time- independent Hamiltonian TL: 



m)) 



-iHt/h 



a|O) + /0|l)) a ® |0> rost ® |0) 6 . (10) 



The quantum information will propagate through the 
network (the "rest"), potentially to emerge in qubit b 
at time T: 

|#(t = T)) « \* f ) = |0) o <g> |0) rcst ® (a\0) + /3|l)) b . 

(11) 

Following the work of Bose the dynamics is of Eq. @ 
is restricted to the states |0) = |0 . . . 0) (the ground state) 
and the N first-excited states |x) = X x \0), where X x is 
the Pauli operator for the qubit at site x. Specifically, 
starting with the initial state |^(0)) = a|0) + (3\a), the 
time evolution given by TL yields the final state 

|*(t))=ae^|0)+/3^/ x (t)| : r) ) (12) 

X 

with 4> — \ Wjt anf i 

f x (t) = e^(ex P H(u; + fi)t]) Xia . (13) 

Here to is a diagonal matrix (with elements ujj ) and Q the 
coupling matrix given by Eq. (|3|). Perfect state transfer, 
as described above, will occur if \fb(T)\ — 1. Note that 
this expression shows that for this type of state transfer 
we need only work with the effective Hamiltonian for this 
subspace, with the matrix form TL e ff = h(uj + Q). We 
now show, using the results of Christandl et al.^ how the 
hypercube network can implement perfect state transfer 
from corner to corner. 

Letting a and b be the (labels for the) qubits at any 
two corners of the hypercube, as in Fig. 1(a), we first set 
LUj = luq in the qubit Hamiltonian ([2]) and approximate 
(J3j> by Cljk ~ (,LUoAjk/2. Then, evaluating f x (t) by using 
the tensor-product form of the adjacency matrix given 
by Eq. (gj) we find 



/*(*) 



= e*V 



-^ t (-z) d [sin(Cc t/2)]' i . (14) 



Ignoring the overall phase <fi, we find that perfect state 
transfer occurs at the transfer time T = tt/^ujq) with 



a = a and j3 = 



-iwoT 



(—i) d p. Note that for phase 



qubits, the initial state can be prepared and the final 
state can be verified using the state tomography tech- 
niques demonstrated by Steffcn et al£Q 

By dynamically switching the qubit frequencies into 
and out of resonance, this state transfer protocol can be 
extended to any pair of qubits in the hypercube. This is 
done by setting ojj = lo^ only if (J, k) is an edge of the 
subcube of the hypercube with a and b as corners, and 
\ujj — Uk\ > C w j,fc otherwise. That is, this network can be 
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"programmed" by tuning the qubit frequencies into reso- 
nance for subcubes of the hypercube. As there are many 
disjoint subcubes,— this network allows for the transfer of 
quantum information in parallel, all with similar transfer 
times. This dynamical switching can be achieved by low- 
frequency control of the phase qubit 's control circuital 
and was demonstrated in a recent experiment.— 

This network can be used for many tasks in quantum 
information processing. By coupling each node of this 
network to an additional set of "memory" qubits, this 
scheme allows for the parallel execution of multi-qubit 
gates between a large (2 d ) number of qubits. A more 
promising application (from the perspective of fault tol- 
erance) builds on this idea to quickly distribute entangle- 
ment between multiple nodes that can then be purified 
locally, and refreshed during computation. This entan- 
glement can then be used for logic gates, error detection, 
or other teleportation-based tasks.— 

IV. HIGHER-ORDER COUPLINGS 

There are many limitations to the success of this 
scheme, which we now analyze. The qubit and rotating 
wave approximations have been used to obtain Eq. ([2]). 
The qubit approximation fails if higher energy states are 
excited. This could occur, for example, during paral- 
lel operation of subcubes, but can be avoided by suffi- 
ciently large detuning.— The rotating wave approxima- 
tion should hold for sufficiently small couplings £ < 0.1^ 

A more significant error is due to the higher-order 
terms in the coupling matrix £1 in Eq. ([3]) . State transfer 
succeeds by engineering the spectrum of the first-excited- 
state subspace. This spectrum is determined by the ad- 
jacency matrix, which for the hypercube can be mapped 
to an angular momentum operator on the vertex space: 
A = 2J X J£ Using this mapping for the phase-qubit hy- 
percube, the effective Hamiltonian H e ff = M w + ^) can 
be written as 

H eff -fi^o(C^ + 2C 2 4 2 + •••), (15) 

where we have ignored an overall constant (7iu>o)- While 
J x has a linear spectrum (leading to perfect transfer) , the 
higher-order terms add quadratic (and higher) terms to 
the spectrum. These terms lead to an intrinsic dephasing 
of the time evolution. 

However, by varying the frequencies of each qubit, 
one can correct for these terms by the following scheme. 
First, label the rows of the d-dimensional hypercube from 
to d (these labels are equal to the Hamming distance 
of the binary label of each node from x a = 0---0). 
Then, setting the frequencies for each qubit in row k to 
u)o + ^o/?2(fc — d/2) 2 leads to the new effective Hamilto- 
nian: 

Heff w hu; {foJ 2 z + (J* + 2( 2 4 + •••). (16) 

One can then show, using perturbation theory in the 
eigenstates of J x , that the choice 02 — 4£ 2 corrects for 




C d 

FIG. 2: Error in state transfer due to higher-order couplings, 
(a) The error 1 — F as a function of coupling strength f, 
where F is the state transfer fidelity for a d = 10 hypercube 
with (circles) and without (squares) the quadratic frequency 
shift. The upper (dashed) line is the function 7r 2 d 2 £ 2 /2, while 
the lower (solid) line is 37r d 3 C 4 /8 + 7r 2 d 6 C 6 /8; these are found 
from a perturbative solution, (b) The error 1 — F as a function 
of hypercube dimension d with £ = 0.005. 



the contribution of J 2 to the spectrum. Thus, by simply 
engineering a quadratic frequency shift across the net- 
work, this lowest order error can be eliminated [sec Fig. 
1(d)]. 

We now demonstrate the effectiveness of this scheme 
by calculating how the error 1—F scales as a function of ( 
and d, shown in Fig. [2] Here F = |(* / |e _i ' H ^' r / fi |* i )| 2 
is the fidelity calculated using the effective Hamiltonian, 
simulating state transfer with the initial state specified 
by a = and (3 = 1. Without any correction, the 
error scales as d 2 ( 2 , shown by the dashed lines, while 
the correction described above yields residual errors pro- 
portional to d 3 ( A and d e ( e , shown by the solid lines. 
These curves are found by calculating the effect of the 
quadratic and cubic terms of the eigenvalue spectrum on 
the state transfer fidelity. These terms arise from both 
the higher-order couplings omitted above (e.g. ( 3 A 3 ) 
and the imposed variation of the qubit frequencies. As 
these qubit frequencies can also be controlled dynami- 
cally, these higher-order errors could be eliminated by 
using well-known quantum control techniques^ Never- 
theless, this single correction reduces the error over two 
orders of magnitude, and becomes tolerable even for large 
(d ~ 10 — 20) hypercube networks. 



V. DECOHERENCE 

Any experimental demonstration will also encounter 
the effects of decoherence and disorder. We model the 
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first by a master equation for the density matrix: 

dtp = -i[H/h, p] + Tr 1 {vjpvj - l{ajaj,p}) 
+\Y, j T; 1 {Z j pZ j -p), 



with af = (Xj ^iYj)/2. Here we have introduced inde- 
pendent energy decay (Ti, also called amplitude damp- 
ing) and dephasing (TU, also called phase damping) pro- 
cesses for each qubit^S For weak decoherence [Ti, ^> 
T = 7r/(£a;o)] we now perturbatively solve this master 
equation given the initial condition p(0) = l^) (vF^. 

We first note that, given the initial condition, p(t) can 
be written as 

pit) = po,o(*)|0)(0|+E x [ft),x(*)|0)(a:H-Px,o(*)|a!)(0|] 

+Ex,„p«.»(*)k)(»l- 

(18) 

Note that p 0) o(0) = H 2 , p ,x(0) = a(3*5 a ^ x , p xfi (0) = 
(x*f3S a ,x, and p x ,x(Q) = \f3\ 2 8 a , x - The density matrix ele- 
ments satisfy the differential equations: 



(17) 



Po,o = Tf 1 J2 X P^,x = Tf 1 (l - p ,o) 

A), a; = ii^xPO,x + * E a ^x,yP0,y — T 2 po,:r 

At,0 = —i^xPxfi — i Ej, ^x,yPy,0 — T 2 Px,0 

Px,y — ^ 5j >,2 Pz.y P:c,z^z,y) 



(19) 



&x,y)- 



-^X ^x.yPx.x 2T 2 /?a;,y(l 

The first equation in (fT9|) can be solved directly to yield 
/9o ,o(«)=l-| ( 3| 2 e-* /Tl . (20) 

Using the conditions lu x — and fl x , y ~ (uoA x , y /2 as 
before, we can use the results of Section III to solve the 
second and third equations in (TT9")) (with x = b) 

po, b (t) = a/?*e-^ e W( +i )d[ sin ( CcJot/ 2)]d 
p 6j0 (t) = a*f3e- t / T2 e~ t « J « t (~i) d [sm(Cu t/2)} d , { ' 

where T^ 1 = (2T 1 y 1 + T7 1 . 

The equation for pi,^, however, cannot be solved ex- 
actly. A perturbative solution can be found, but requires 
a careful examination of the eigenstates of Q and how 
they are mixed by the decoherence terms. After a long 
calculation-^ we find 



p b , h {t) = |/3| 2 e-*/ Tl ^(-l) d -"g„(t)cos(C^(d-n)), 

(22) 



n=0 



where we have defined the functions 

_ d!(2-<5» d )2"- 2d - 2 P _ A t 
h PKn-2p)!(d-n+p)! e 



and the decay constants 
2 



-Lrh 



(d-n + 2p)\ 
p\(d — n + p)\ 



(23) 



(24) 
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FIG. 3: Probability pb,b{t) as a function of time for state 
transfer with /3 — 1. The symbols correspond to numerical 
simulations of the master equation with u>o/2tt = 5 GHz, 
£ = 0.005, T\ — — 120 ns, for hypercubes of dimension 
d — 2 (circles) and d = 10 (squares). The corresponding lines 
are evaluations of the perturbative solution given by Eq. (I22|) . 



As shown in Fig. [3J this solution for pb,b(t) compares well 
to direct numerical simulations of the master equation 
(|17p as well as to more detailed simulations (not shown) 
of a master equation for the phase qubit system involving 
a large number of non-qubit states (e.g. hundreds of 
states for d = 3). 

Using these results we can calculate the state trans- 
fer fidelity F = (*/|p(T)|*/), at the transfer time 



T = tt/(Cw ) with |*/) = a|0) + (-i) d e 



-iu T 



P\b): 



F 



<|4 e -T/Ti 



-2|a| 2 | 

E„=o#»( T )- 



-T/T 3 



(25) 



An important feature of this solution is that, so long as 
the frequency loq and the coupling £ remain fixed, the 
state transfer fidelity has a lower bound that is indepen- 
dent of the size of the network. To show this, we follow 
Bos&i 2 - and consider the fidelity F avg averaged over the 
Bloch sphere (the worst case can be bounded similarly): 

F avg = \ V T/Tl + \e- T ^- + \e~ T ^ £>(T). 

' ' ' n=0 

(26) 

Using the fact that X pn < 2/T^, we find that 
J2Lo9n(T)>e- 2T / T * and thus 



1 1 

F > p 

avg ~ 2 6 



-T/Ti 



-T/Ti 



-2T/T 2 



(27) 



Since the transfer time T = %/ (C^o) is independent of d, 
so is this bound on the fidelity. Furthermore, this simple 
formula shows that the fidelity can be greater than the 
classical limi^ of 2/3 for existing experimental devices. 
For example, with ujq/2-k — 5 GHz, Q — 0.005, and the 
decoherence times T\ = 120 ns and T2 = 80 ns,— we find 



T 



20 ns and thus F avg > 0. 
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0.2 0.4 0.6 0.8 1 

FIG. 4: Fidelity of state transfer due to disordered couplings. 
Each point for d — 6 (circles), d — 10 (squares), and d — 16 
(triangles) represents the fidelity averaged over many (10000, 
1000, and 100, respectively) realizations of disordered cou- 
plings (see text). The corresponding lines (solid, dashed, and 
dotted) are fitted Gaussians. 



Note that this bound on the fidelity is specific to our 
decoherence model. If the qubit decoherence times de- 
pend on the size of the circuit, this model should be 
improved. Physically, the complexity of the wiring de- 
sign (the cross-overs and interconnects) could introduce 
significant loss using traditional dielectrics^ However, 
advanced fabrication using vacuum gaps should remove 
this potential decoherence source^ Note also that for 
fixed capacitors (instead of fixed uiq and £), there is 
in fact a weak dependence of the transfer time on d: 
T(d) fa T(l)[l + (d- 1)C] 3/2 - However, for the exper- 
imental parameters given above, this dependence can be 
neglected for d < 20. 

VI. DISORDERED COUPLINGS 

No experimental circuit will have identical qubits — 
there will always be some variation. This variation in 
parameters is equivalent to disorder, and in this context 
we might expect a blocking of transport due to Anderson 
localization. 34 That is, by mapping propagation through 
the hypercube to an effective spin, the information trav- 
els along a one-dimensional path (row by row)J£ This 
was studied in previous work^ in which static disorder 



was introduced into the one-dimensional effective Hamil- 
tonian for perfect state transfer. There it was found that 
disordered couplings play a much larger role than disor- 
dered qubit frequencies, with results in agreement with 
localization theory. Note that we have already allowed 
the qubit frequencies to be under experimental control, 
so here we only consider disordered couplings. 

To simulate the effects of static disorder on multi- 
dimensional state transfer, we use use the full subspace 
Hamiltonian and introduce disorder into the coupling 
matrix by Qjk = wq(C + z jk)Ajk/2 (the higher-order 
terms are neglected here) where each Zjk is a random vari- 
able uniformly distributed over the range — A£ < Zjk < 
A^. The transfer fidelity (with a = and (3 = 1) is then 
calculated and averaged over this uniform distribution. 
The averaged fidelity is shown in Fig. These results 
are consistent with F ~ e~ d / e , where i ~ c(A()~ 2 . This 
implies that, for very large d, transport is exponentially 
suppressed due to localization. However, even a hyper- 
cube of modest dimension, where localization effects are 
small, can accomodate a large (2 d ) number of qubits. For 
example, while state-of-the-art fabrication can achieve 
variation of order 1% or less, Fig. Q] shows that even 
for 10% variation in a network with d = 10 one still has 
F > 0.95. 



VII. CONCLUSION 

In conclusion, we have proposed that the hypercube 
network can be used to route quantum information be- 
tween superconducting qubits. Using tunable phase 
qubits, this network can be programmed to achieve high- 
fidelity quantum state transfer between any two nodes 
and allows for rapid distribution of entanglement in a 
superconducting quantum computer. The dominant er- 
ror mechanisms have been analyzed, including higher- 
order couplings, decoherence, and disorder. Our results 
indicate that using existing qubit fabrication and experi- 
mental methods, a demonstration of these operations can 
be performed. More generally, these initial results moti- 
vate continued study of novel quantum transport in these 
artificial solids. 
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